Numerical study of finite size effects in the one-dimensional two-impurity Anderson 

model 
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We study the two-impurity Anderson model on finite chains using numerical techniques. We 
discuss the departure of magnetic correlations as a function of the interimpurity distance from a 
pure 2fc_F oscillation due to open boundary conditions. We observe qualitatively different behaviors 
in the interimpurity spin correlations and in transport properties at different values of the impurity 
couplings. We relate these different behaviors to a change in the relative dominance between the 
Kondo effect and the Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction. We also observe that 
when RKKY dominates there is a definite relation between interimpurity magnetic correlations and 
transport properties. In this case, there is a recovery of 2fc^ periodicity when the on-site Coulomb 
repulsion on the chain is increased at quarter-filling. The present results could be relevant for 
electronic nanodevices implementing a non-local control between two quantum dots that could be 
located at variable distance along a wire. 

PACS numbers: 72.15.Qm, 73.23.-b, 73.63.Kv, 75.40.Mg 
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I. INTRODUCTION 

It is well-known that the main properties of the two- 
impurity Anderson model (TIAM)- are determined by 
the competition between the Kondo effect^ generated by 
each impurity and the Ruderman-Kittel-Kasuya-Yoshida 
(RKKY) interaction between the impurities Renewed 
interest in this model comes from the observation of the 
Kondo effect in electronic devices where a nanoscopic rel- 
atively isolated region, a quantum dot (QD), can act 
as a single spin-l/2£ In particular, it has been ob- 
served in a nanoelectronic device formed by two QDs 
coupled through an open conducting region^ that trans- 
port through one of the QDs depends on the occupancy 
and coupling of the other QD via the RKKY interac- 
tion. The main issues which are involved in this and re- 
lated experiments 7 - are the possibility of tuning different 
regimes dominated by the Kondo effect or by the RKKY 
interaction, and the dependence of transport properties 
with the sign of the magnetic correlation between the 
impurities. 

Analytical studies of the TIAM are usually concerned 
with long-distance behaviors^ However, it should be 
emphasized that the properties of mesoscopic devices 
are determined by finite size effects rather than by bulk 
physics. Finite size effects have been addressed in a study 
of a double QD in an Aharonov-Bohm ring£ and it was 
found that RKKY dominates the transport properties. 
More closely related to the problems we consider in this 
work is a study of two impurities coupled by a finite 
one-dimensional (ID) wire^ where it was shown that the 
RKKY interaction is always dominant due to the strong 
reduction of the Kondo temperature by finite size effects. 
In another slave-boson studyi^ it was concluded that the 
presence of different transport regimes depends on the 
sign of the RKKY interaction. In these two studies how- 
ever, the impurities are located at the edges of the chain 
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FIG. 1: (Color online) Picture of model (JU, L = 12, R = 3. 



and hence the effect of the outer sites in Fig.[TJ or "leads" , 
were not considered. 

The main motivation for the present study is then to 
understand finite size effects in the ID TIAM which can 
be considered as a model of a device formed by two QDs 
coupled through a conducting chain and connected to ID 
leads (Fig.[T]). Another motivation for this work is to un- 
derstand the effects of electronic correlations which are 
present in materials such as carbon nanotubesii used in 
experimental realizations of quasi-lD devices^ Coulomb 
correlations in the chain affect both the Kondo temper- 
ature and RKKY interaction^^ 

In this work we show numerically that the spin-spin 
correlation between the two impurities as a function of 
the interimpurity distance R, S(R), for a fixed length 
L of the whole chain, departs from the predicted 2kp 
long-distance behavior for the TIAMi&i We show how 
a regime where the RKKY interaction dominates over 
the Kondo effect appears as the impurity couplings are 
decreased. We also provide a measure of the conductance 
G which is more relevant for applications to nanodevices. 
We show that this quantity, particularly in the regimes 
where RKKY dominates, follows an oscillatory behavior 
with R at a fixed L related to that of S(R). We also 
show that an on-site Coulomb repulsion on the chain, at 
quarter-filling, restores the 2fcp periodicity of S(R) and 
G(R). Finally, we provide indications that most of these 
finite size effects found are due to the open boundary 
conditions adopted for our system. 



The paper is organized as follows. In Sec. |TT] we de- 
fine the Hamiltonian model and we describe the method 
employed in the calculations. In Sec. Mil we present re- 
sults for S(R) and G(R), first for noninteracting leads 
and then in the presence of correlations on the chain. 
We conclude in Section IIVI by providing a summary of 
the results and by suggesting possible connections with 
other theoretical approaches as well as with experiments. 

II. MODEL AND METHODS 

The ID two-impurity Anderson model with Hubbard 
repulsion on the chain is defined by the Hamiltonian: 

t ^ ] ( c r 1 +fj,a- c ri<J + c r 2 +/jff C i'20 + H.C.) 
fl — ±l,cr 

+ U' ^ n^n^+e' (n lT +n a ) (1) 

i— ri,r2 i-ri,r2 

where the notation is standard. The Anderson impurities 
or QDs, with parameters U' , e', are symmetrically located 
with respect to the center of the chain, i.e. at sites r 12 = 
(L±R + l)/2, with R odd (see Fig. Q]). The QDs' are 
connected to the rest of the system A with a hopping t'. 
The subsystem A comprises the leads (i < n, i > r^) and 
the region between both impurities (n < i < 1*2), and it 
is described by a Hubbard Hamiltonian with couplings to 
and U . to = 1 is adopted as unit of energy. 

Model {1} will be studied on finite clusters using 
density matrix-renormalization group (DMRG), 1 - Open 
boundary conditions (OBC) were adopted throughout. 
We would like to emphasize that OBC are the realistic 
boundary conditions for the devices that motivate the 
present work. There has been a previous study of a re- 
lated model using DMRG but the two impurities were 
fixed at the chain edges^ that is without leads. 

We provide a measure of the response of the system to 
the application of a small bias voltage AV = Vr — Vl 
(Vr = —Vl), where Vl (Vr) are on-site potentials ap- 
plied to the ten sites at the edge of the left (right) lead. 
Although these ID leads have to be connected to mas- 
sive two- or three-dimensional metallic or semiconducting 
contacts, we believe that as a first step it is necessary to 
consider the system formed only by the two QDs con- 
nected to the ID leads. The current J hit) (Jn(t)) on the 
bond connecting the left (right) QD to the left (right) 
lead is computed with the time-dependent DMRG^ In 
the case of a single QD and non-interacting leads, this 
numerical setup has been show n 17 ' 18 i 19 ' 20 to reproduce 
the results for equivalent systems which were treated 
analytically. 21 In the following, we employed the "static" 
algorithmic which although less precise than the "adap- 
tive" algorithm 18 still gives qualitatively correct results, 
as we have explicitly checked in some cases ; 22 ' 23 but is 
much faster than the adaptive scheme thus allowing to 
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FIG. 2: (Color online) (a) Spin-spin correlations between the 
impurities, and (b) conductance, as a function of the interim- 
purity distance R, at n = 1, for t' — y/2/2 (circles) and t' = 1 
(squares), on the L = 120 chain except otherwise stated. In 
(a), dashed lines show fits to 2fcp oscillations. 

explore different couplings and densities and a wide range 
of R and L. All the results reported below correspond 
to AV = 0.01. We compute J(t) = (J L {t) + J R (t))/2, 
which follows with time very approximately a sinusoidal 
function, and we adopt as a measure of the response or 
"conductance" , G, the average of the sinusoidal fitting to 
J(t)/AV over a half-period. This qualitative behavior is 
mostly independent of various criterions that have been 
proposed to quantify G*£ We want to emphasize that 
our main purpose in this study is to obtain relative val- 
ues of G as a function of the interimpurity distance R. In 
this sense, the application of the bias potential on only 
ten sites at the edges of the system instead of applying 
to the whole leads is then a necessary condition to treat 
all distances R on an equal foot. 

We adopted V = 8 and e' = —4 (symmetric point) in 
order to ensure a strong magnetic character of the impu- 
rities which we checked in all cases. In the following we 
will consider two values of t 1 , \/2/2 and 1. The Kondo 
effect for the single impurity Anderson model for t' = 1 
was studied using the same method as in the present 
work in Ref. [22|. Most of our calculations reported be- 
low were obtained for L = 120. The truncation error in 
DMRG calculations for this cluster was kept below 10~ 6 . 
In all cases we worked in the subspace of total S z = 0. 
Thus, taking into account the isotropy of the Hamilto- 
nian (P) we computed the spin-spin correlations between 
sites i ^ j as Sij — 3(5f In particular, the spin-spin 
correlation between impurities is S(R) — 5(S^ 1+R S Z 1 }. 



III. RESULTS AND DISCUSSION 

Let us start by discussing results at half-filling, n = 1. 
In Fig. [Ha), it can be seen that for t' = \/2/2 and t' = 1, 
S(R) presents a 4-site period which corresponds to a kp 
oscillation. To be more precise there is a kp modula- 
tion on top of the antiferromagnetic (AF) 2kp oscillation. 
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FIG. 3: (Color online) (a) Spin-spin correlations between the 
impurities, and (b) conductance, as a function of the interim- 
purity distance R, for t' = V2/2, at n = 0.5, parameterized 
with U, on the L = 120 chain except otherwise stated. 

This 2kp component can be fitted for R odd by the law 
-1.102/a; - 59 0' = x/2/2) and -1.021/x 1 38 (*' = 1). The 
kp component decays much slower, as x~ 0,22 (t' = y/2/2) 
and £-°- 68 (f = 1). At this point we could advance the 
hypothesis that the departure of the present results from 
the pure 2kp oscillation is due to the OBC used in DMRG 
calculations. That is, 2kp oscillations starting from the 
chain edges would modulate the magnetic correlations 
from each of the impurities and hence, as a function of 
R, this would enhance the kp component of S(R). The 
large intensity of S(R) and its survival at long distances 
for t' = V2/2 suggest that we are in the regime where 
the RKKY interaction dominates over the Kondo effect. 
On the other hand, the smaller amplitude and rapid sup- 
pression of S(R) with R for t' = 1 indicate the presence 
of a strong screening of the magnetic moment of the im- 
purities implying a dominance of the Kondo effect. 

The conductance (Fig. Efb)) presents an oscillation 
with R that follows that of S(R). There is a set of min- 
imum values of G(R) for R = Am + 3 (m integer) that 
correspond to the stronger AF spin-spin correlations be- 
tween impurities, and a set of maximum values of G(R) 
which occur at R = Am + 1 corresponding to the weaker 
AF S(R). The expected difference in amplitude of G(R) 
between both values of t' considered, since the current is 
proportional to t', can also be observed. The most im- 
portant feature in these results is that, for the case of 
t' = V2/2, there is a factor of 2 or larger between the 
conductance for R = Am + 1 with respect to the one for 
R = Am + 3. We believe that this large difference can 
be detected experimentally in appropriate devices. In 
our calculation, the maximum of G(R) could not possi- 
ble reach the expected unitary limit as a consequence of 
applying the bias potential to only few sites on the leads 
as discussed in Section [TT1 

Although half-filling is perhaps more realistic for ac- 
tual devices, it is interesting to study also the case of 
quarter-filling, n = 0.5. In this case it is possible to study 
transport properties in the presence of electron correla- 



FIG. 4: (Color online) (a) Spin-spin correlations between the 
impurities, and (b) conductance, as a function of the interim- 
purity distance R, for t' = 1, at n = 0.5, parameterized with 
U, on the L = 120 chain. In (b) the singlet-triplet energy 
Asr/t' is shown for t' = 1 (stars) and t' — y2/2 (crosses). 

tions on A. In Fig.[3^a) , for t' = V2/2, it can be seen that 
for the noninteracting chain (U = 0) S(R) presents an 
oscillation with period 8 which again corresponds to an 
overall kp oscillation. As in Fig. HJa), results for L = 80 
show an oscillation in S(R) slightly larger in amplitude 
to that for L = 120. A systematic study of the effect of 
L is presented at the end of this work. Similarly to what 
happened at n = 1, the conductance, shown in Fig. E^b), 
has the same dependence with R as the spin-spin correla- 
tions. At n = 0.5, the minimum (maximum) of G occurs 
at R = 8m + 5 (R = 8m + 3) which coincides with the 
values of R at which S(R) has the strongest (weakest) 
AF values. 

Let us turn on the interaction on A. It can be seen 
in Fig. [3]Ja) that there is a rapid suppression of the fcp 
modulation in S(R) with increasing U leaving behind a 
well-defined 2kp oscillation which is already apparent at 
U = 2. For U > 2 the spin-spin correlations at R — Am + 
3 become slightly positive, that is ferromagnetic (FM). 
Consequently, as it can be observed in Fig.[3Jb), there is 
a similar change in the dependence of the conductance 
with R, with the minimums (maximums) of G located at 
the same values as the AF (FM) correlations in Fig.^a). 
Notice that for U > 2 the ratio between the maximum 
and the minimum values of G becomes equal to 8 or 
larger. Again this very large ratio should be detectable 
experimentally if a device with the geometry of Fig. [1] 
and working at quarter-filling could be fabricated. 

In Fig. HJa), we show the spin-spin correlations be- 
tween impurities at n = 0.5 for t' = 1. At this filling, 
there is an additional important difference with the cor- 
responding results for t' = y/2/2 depicted in Fig. [3Ja) 
which is that the 2kp periodicity is now clearly domi- 
nant. This behavior of S(R) supports the idea that these 
two values of t' belong to different regimes. On the other 
hand, the behavior of the conductance with R shown in 
Fig. BJb), m the noninteracting case U = 0, is similar to 
the one for t' = \2/2, shown in Fig. [3Jb), with a large 



kp component. However, in the presence of interactions 
on the chain the amplitude of G(R) is suppressed by in- 
creasing U > for t' = 1, while it seems enhanced by 
U for H = y/2/2. This suppression of G(R) would be 
expected in a Luttinger liquid since it is proportional to 
K p which decreases with increasing U i 22 i 24 For t 1 = 1 , 
it can be also observed in Fig. QJb) that only for the 
largest value of the interaction considered, U = 4, a 2kp 
periodicity becomes dominant. Finally, the relation be- 
tween the location of the maximum and minimum values 
of G(R) with the ones of S(R) observed in Fig. ^[b) is 
absent for small U but it is recovered for U = A. 

These clearly different behaviors of S(R) and G(R) 
for the two types of impurities considered have to be 
traced to the essential competition between the Kondo 
effect and RKKY interactions. For the noninteracting 
case, the Kondo temperature may be estimated as 2 Tk — 
t'y/U'/(2t ) exp {-irt U'/{8t' 2 )) (in units of k B = 1) 
which gives T K = 0.0026 for t' = y/2/2 and T K = 0.086 
for t' — 1. That is, if for the case of t 1 — \[2j2 it is rea- 
sonable to assume that RKKY dominates over the Kondo 
effect, since Tk for t' = 1 is 33 times larger while the 
effective RKKY interaction^ K. ~ t' 2 would be only a fac- 
tor of 2 larger, then it is plausible that the case of t' = 1 
belongs to the regime where the Kondo effect dominates 
over the RKKY interaction. Particularly for the case of 
f = V2/2 where RKKY can be assumed as dominant, it 
is instructive to estimate JC as the difference in ground 
state energy between the singlet and triplet states. For 
this case one gets that Ast oscillates with R with max- 
imum values of K, ~ Ast ~ 0.04 (n = 1). This value is 
larger than the above estimated Tk for this impurity. At 
n = 0.5 the oscillation of Ast with R has period 8 as it 
can be observed in Fig. B|b). Although these estimates 
are instructive and consistent with our numerical results, 
one should bear in mind that the above expression for Tk 
is valid in the bulk limit while the results for Ast were 
obtained for the L — 120 chain. The strong suppression 
of Tk due to the finite size of the system found in Ref. 9 
may not appear in our system due to the presence of the 
leads which may provide room for the Kondo cloud to 
develop. We would also like to notice that in all cases we 
examined the ground state is a singlet, which is consis- 
tent with the result that S(R) is always AF or at most 
weakly FM. 

Finally in Fig. \5\ we show the dependence of S(R) at 
a fixed interimpurity distance R = 5 with the length 
L. The period of S L (R = 5) is 4 (n = 1) and 8 (n = 
0.5) for the two types of impurities considered, which 
correspond to wave vector kp. It is interesting to examine 
here the dependence of these results with the boundary 
conditions (BC) of the system. Periodic BC (PBC) and 
antiperiodic BC (APBC) have, separately, a period of 8 
as with OBC. However, if for each L we adopt the BC 
with minimum energy, then a doubling of the wave vector 
of the oscillation is obtained. Using exact diagonalization 
we obtained that the energy is minimum for PBC for 
L = 4m (m integer) and for APBC when L = 4m + 2. 
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FIG. 5: (Color online) Spin-spin correlations between the im- 
purities at 7? = 5, as a function of the total length L, for 
t' = \/2/2 (circles), t' = 1 (squares), and (a) n = 1, (b) 
n — 0.5. The inset shows the scaling of Sl(R) with L, for 
t' = V2/2, n = 0.5, and for various values of R as indicated. 

The minimum of energy for L = Am + 1 and L = Am + 3 
falls alternatively on PBC or APBC depending on R. 
The resulting Sp(R = 5) is shown in Fig. OJb). This 
result is then another example of the effects caused by 
OBC. 2 - The inset of Fig. [5] shows a slow dependence of 
Sl(R) with L for various R. Since, particularly for t' = 
y/2/2, G(R) follows the oscillation of S(R), it is quite 
likely then that an oscillation of the conductance could 
be experimentally observed in a real device for a fixed 
interdot distance as the length of the ID leads is varied. 



IV. CONCLUSION 

In conclusion, we have shown for the two-impurity 
Anderson model on finite chains, the presence of differ- 
ent behaviors of the interimpurity magnetic correlations 
S(R) for the two values of the coupling between the impu- 
rities and the chain, t 1 . We suggested that these different 
behaviors indicate that the relative dominance between 
the Kondo effect and the RKKY interaction can be tuned 
by the single parameter t' . In the case when RKKY dom- 
inates we found important oscillations in the conductance 
with interimpurity distance, G(R), following that of the 
spin-spin correlations. It would be tempting to relate 
these different behaviors of S(R) and G(R) to the pres- 
ence of an unstable fixed point suggested by numerical 
renormalization group calculations. 3 Although certainly 
such fixed point could not strictly appear in a finite-size 
calculation some traces of its presence, for example an 
admixture of the singlet and triplet states, could be de- 
tected. However, we found that the ground state of our 
model is always a singlet consistently with dominating 
AF interimpurity correlations. It would be interesting 
to explore other parameters of our model to find such a 
critical points 

Finally, we have observed that the kp modulation of 
the 2kp oscillation of S(R) and G(R) observed for non- 
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